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Abstract 

The density matrix renormalization group method is used to examine possi- 
bilities of electronic ferroelectricity in the spinless Falicov-Kimball model. The 
model is studied for a wide range of parameters including weak and strong 
interactions as well as the symmetric and unsymmetric case. In all examined 
cases the < d + f >-expectation value vanishes for vanishing hybridization V, 
indicating that the spinless Falicov-Kimball model does not allow for a ferro- 
electric ground state with a spontaneous polarization. 

PACS nrs.: 71.27.+a, 71.28.+d, 71.30.+h 
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The Falicov-Kimball model (FKM) has become, since its introduction JIJ in 1969, 
one of the most popular examples of a system of interacting electrons with short- 
range interactions. It has been used in the literature to study a great variety of 
many-body effects in metals, of which metal-insulator transitions, mixed-valence 
phenomena, and charge- density waves are the most common examples B|. In the 
last few years the FKM is extensively studied in connection with the exciting idea 
of electronic ferroelectricity || §], |5], || f?|. It is generally supposed that the ferro- 
electricity in mixed- valent compounds is of purely electronic origin, i.e., it results 
from an electronic phase transition, in contrast to the conventional displacive fer- 
roelectricity due to a lattice distortion. Since the FKM is probably the simplest 
model of electronic phase transitions in rare-earth and transition-metal compounds 
it was natural to start to examine possibilities of electronic ferroelectricity just on 
this model. An important milestone on this way was the recent work by Portengen at 
al. ||, Q. They studied the FKM with a k-dependent hybridization in Hartree-Fock 
approximation and found, in particular, that the Coulomb interaction U between the 
itinerant (^-electrons and the localized /-electrons gives rise a non-vanishing excitonic 
< d + f >-expectation value even in the limit of vanishing hybridization V —> 0. As 
an applied (optical) electrical field provides for excitations between d- and f-states 
and thus for a polarization expectation value Pdf =< dffi >, the finding of a spon- 
taneous Pdf (without hybridization or electric field) has been interpreted as evidence 
for electronic ferroelectricity. This result motivated further theoretical studies of the 
model. Analytical calculations within well controlled approximation (for U small) 
performed by Czycholl || in infinite dimensions do not confirm the existence of elec- 
tronic ferroelectricity. In contrast to results obtained by Portengen et al. [0, |4j] he 
found that the FKM in the symmetric case (Ef = 0,nf = = 0.5) does not allow 
for a ferroelectric ground state with a spontaneous polarization, i.e., there is no non- 
vanishing < d + f >-expectation value in the limit of vanishing hybridization. The 
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same conclusion has been also obtained independently by extrapolation of small- 
cluster exact-diagonalization calculations in the one dimension for both intermediate 
and strong interactions M. In these regions the finite-size effects are negligible and 
thus the results can be satisfactory extrapolated to the thermodynamic limit. Very 
recently Zlatic et. al. J7| analysed the static susceptibility for spontaneous polariza- 
tion in the ordinary FKM (V = 0) using the exact solution of the model in infinite 
dimensions. They found, in particular, that the spontaneous hybridization suscep- 
tibility diverges at T — > for U = and probably for nonzero U too, indicating a 
possibility of a nonzero < d + f >-expectation value in the ground state. 

In order to shed some light on these controversial results we performe an exhaus- 
tive study of the model for a wide range of model parameters. While in our previous 
paper the possibilities for the electronic ferroelectricity have been studied only in the 
symmetric case and intermediate and strong interactions || we extend here calcula- 
tions also to the unsymmetric case and small interactions. To suppress the finite-size 
effects in the weak-coupling limit one has to consider much larger clusters than can 
be treat numerically within the exact-diagonalization calculations [L ~ 12). For 
this reason we have decided to use in numerical computations the density matrix 
renormalization group (DMRG) method || that allows to treat several times larger 
clusters and still to keep the high accuracy of computations. 

As mentioned above the FKM is based on the coexistence of two different types 
of electronic states in a given material: localized, highly correlated ionic-like states 
and extended, uncorrelated, Bloch-like states. It is accepted that insulator-metal 
transitions result from a change in the occupation numbers of these electronic states, 
which remain themselves basically unchanged in their character. Taking into ac- 
count only the intra-atomic Coulomb interaction between the two types of states, 
the Hamiltonian of the spinless FKM with hybridization can be written as the sum 
of four terms: 
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H = £ Ujd+dj + ftMdi + E f ]T + V E + h-c, (1) 

ij i i i 

where /+, /j are the creation and annihilation operators for an electron in the local- 
ized state at lattice site % with binding energy Ef and df, d{ are the creation and 
annihilation operators of the itinerant spinless electrons in the <i-band Wannier state 
at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant d electrons between sites i and j. These intersite hopping 
transitions are described by the matrix elements Uj, which are —t if i and j are the 
nearest neighbors and zero otherwise (in the following all parameters are measured 
in units of t). The second term represents the on-site Coulomb interaction between 
the ci-band electrons with density rid = \ J2i dfdi and the localized / electrons with 
density rif = jj^iftfi^ where L is the number of lattice sites. The third term stands 
for the localized / electrons whose sharp energy level is Ef. The last term represents 
the hybridization between the itinerant and localized states. 

In order to test the convergence of the DMRG method for the FKM with hy- 
bridization let us first consider the case U — (Ef = 0). In this case the Hamiltonian 
(1) can be diagonalized exactly yielding the following expressions for the ground-state 
energy and the expectation value 



E a = ±Ue*-(4 + W*) 1 < 2 ), (2) 

k 



Pdf (e2 + 4V)V2' (3) 

where 6k = \ Y^ij t i je' lk ^ Ri ~ Rj \ These expressions are valid for arbitrary L and thus 
they can be used directly for comparison with DMRG results obtained on finite 
clusters. Such a comparison shows that the DMRG results reproduce the exact 
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results with high accuracy (~ 1CT 5 ) even if a relatively small number of states is 
kept (m = 128), thereby the clusters up to 100 sites become accessible for numerical 
studies. This is illustrated in Fig. la where the difference AE = E e g xact — E^ M 
is plotted as a function of L for m = 128 and two values of V. It is seen that the 
difference AE is smaller than 10~ 5 for all clusters of size L < 100. One can expect 
that the finite-size effects will be considerably reduced for such cluster sizes and thus 
the DMRG results could be satisfactory extrapolated to the thermodynamic limit 
also for weak interactions. Numerical results obtained for the ground-state energy 
E g (Fig. lb) and the Pjf expectation value (Fig. 2a) fully confirm this conjecture. 
It is seen that the finite-size effects are small for both E g and Pdf and so the most 
important feature, and namely the behavior of Pdf for V — > can be reliably deduced 
from these results. In all cases the < d + f >-expectation value vanishes in the limit 
V — * 0, so there is no spontaneous polarization in the one-dimensional spinless FKM 
with nearest-neighbor hopping |9j for U — 0. Using the expression (3) one can show 
immediately that the same behavior holds in any dimension. In Fig. 2b we have 
displayed numerical results for P^ obtained in the thermodynamic limit for d = 1, 2 
and d = oo (Bethe lattice). In all cases Pdf vanishes continuously for V —> 0, which 
seems to contradict recent results obtained by Zlatic et. al. f?[] for the ordinary FKM 
(V = 0) in infinite dimensions. Indeed, a comparison of our d = oo results for the 
< d + f >-expectation value with d = oo results obtained by Zlatic et. al. |7j] for the 
spontaneous hybridization susceptibility Xhyb leads to very interesting observation. 
While the hybridization susceptibility Xhyb diverges for V = 0, the Pdf expectation 
value goes to zero in the limit of vanishing V. This can be shown also analytically 
using the explicit expression (Eq. 3) for the Pdf expectation value. Replacing the 
summation in Eq. 3 by integration over the energy and taking the constant density 
of states, for simplicity, one can easily show that the Pdf expectation value (in the 
limit of vanishing V) behaves like VlnV . Thus Pdf vanishes at V = 0, although Xhyb 
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diverges at this point. This shows that one has to be very careful in doing conclusions 
about the Pdf behavior from the spontaneous hybridization susceptibility. 

The same behavior of the model we expect also for nonzero U. To verify this we 
perform an exhaustive study of the model for a wide range of the model parameters. 
In addition to the half-filled band case n = 1 we examine here in detail also the case 
n = 0.5, Ef < 0, where Zlatic et. al. [0] found a numerical evidence for a spontaneous 
polarization (the divergence of the spontaneous hybridization susceptibility). 

Let us first discuss the DMRG results obtained for in the half-filled band case 
and small but finite U. This case was not discussed in our previous paper || since 
small-cluster exact-diagonalization calculations were not able to give reliable answers 
due to large finite-size effects in this region. As mentioned above within the DMRG 
method much large clusters can be treated numerically, thereby the finite-size effects 
are considerably reduced. The weak-coupling numerical results for Pdf are displayed 
in Fig. 3 for U = 0.5 and several values of Ef. Obviously the one-dimensional FKM 
does not exhibit a ferroelectric ground state with a spontaneous polarization in the 
weak-coupling limit. For all examined values of Ef the < d + f >-expectation value 
vanishes for V — > 0, and it is demonstrated that this result is independent of L. 

The same calculations are performed also for n = 0.5 and Ef < 0. Since the 
numerical results obtained by Zlatic et. al. |7j show on the possibility of a spon- 
taneous polarization in this region a very detailed analysis of the model is done in 
this parameter space for both weak and strong interactions. The typical behaviors 
obtained by numerical simulations are presented in Fig. 4. Again there is no sponta- 
neous polarization, neither for weak nor for strong interactions. In all cases Pdf — > 
for V — ► 0. The finite-size effects are small for both weak and strong interactions 
and thus these results can be satisfactory extrapolated to the thermodynamic limit. 

Thus we can summarize our results. The numerical studies of the spinless FKM 
by density matrix renormalization group method did not confirm a possibility of 
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a ferroelectric ground state with a spontaneous polarization in this model. In all 
examined cases (weak and strong interactions, the symmetric and unsymmetric case) 
the Pdf expectation value vanishes for V — > 0, in contrast to what is expected from the 
behavior of the spontaneous hybridization susceptibility J7J. This fact considerably 
reduces (although does not exclude definitively) the possibility for the existence of a 
ferroelectric ground state with a spontaneous polarization in the spinless FKM. 

This work was supported by the Slovak Grant Agency for Science under grant 
No. 2/7021/20. 
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Figure Caption 



Fig. 1. (a) The difference AE = Ef act - E° MRG as a function of L calculated 
for U = 0,m = 128 and two different values of V. (b) The ground-state energy as a 
function of hybridization for U = and m = 128. Numerical curves corresponding 
to L — 12, 24 and 48 are calculated by DMRG method, while the exact analytical 
solution in the thermodynamic limit L = oo is obtained from Eq. 2. 

Fig. 2. (a) Hybridization dependence of the < d + f >-expectation value calculated 
for finite clusters (DMRG results) and in the thermodynamic limit (from Eq. 3) for 
U = 0. (b) Hybridization dependence of the < d + f >-expectation value calculated 
in the thermodynamic limit (from Eq. 3) for d = 1, 2 and d = oo. 

Fig. 3. Hybridization dependence of the < d + f >-expectation value calculated 
for U = 0.5 and several values of Ef (DMRG results) at half-filling (n = 1). To 
reveal the finite-size effects the calculations have been performed for several finite 
clusters. 

Fig. 4. Hybridization dependence of the < d + f >-expectation value calculated for 
U = 0.5 (a) and U = 5 (b) and several values of Ef (DMRG results) at quarter-filling 
(n = 1/2). 
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